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Abstract 

A  two  dimensional  model  was  constructed  for  a  high  temperature  PBI/H3PO4  fuel  cell,  the  model  was  focused  on  the  cathode  side  since 
it  is  the  most  important  component  of  the  overall  polarization  loss  in  unit  cell.  Linear  sweep  voltammetry  (LSY)  technique  was  used  to  sup¬ 
ply  the  model  with  the  most  important  kinetic  parameter — cathodic  exchange  current  density.  The  cathode  polarization  was  calculated  by  the 
resulting  two  dimensional  model  and  the  ohmic  polarization  was  estimated  by  electrochemical  impedance  spectroscopy  (EIS),  the  overall  polar¬ 
ization  curve  was  simulated  according  to  the  above  procedure;  three  unit  cells  with  different  catalysts  were  assembled  and  their  steady  state 
polarization  curves  were  simulated  based  on  the  mathematical  model.  Good  agreement  between  the  simulated  and  experimental  results  was 
achieved. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Great  success  has  been  achieved  with  proton  exchange  fuel 
cells.  Advanced  features  of  this  technology  include:  higher 
power  density,  environment  friendly  operation,  all-solid  con¬ 
struction  and  therefore  reduced  materials  corrosion,  but  still 
some  serious  challenges,  such  as  the  catalytic  activity  when  the 
operation  temperature  is  below  80  °C;  catalyst  tend  to  be  poi¬ 
soned  by  impurities  (CO,  nitrides  or  sulfides)  in  the  fuel  stream; 
a  major  part  of  the  water  product  is  liquid  and  this  can  result  in 
flooding  of  the  electrode;  still,  if  operated  below  80  °C,  the  small 
difference  between  operational  temperature  and  ambient  tem¬ 
perature  requires  intensive  cooling,  especially  when  the  ambient 
temperature  is  high  or  the  cell  is  under  high  load.  Recently,  many 
researchers  proposed  different  ways  to  solve  the  above  problems, 
their  efforts  include  developing  high  activity  electro-catalysts  to 
increase  cell  performance  [1,2],  improving  the  electrode  struc¬ 
ture  for  higher  catalyst  utilization  [3],  developing  CO  tolerant 
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catalyst  to  alleviate  catalyst  poisonous  effects  [4]  and  proposing 
novel  gas  diffusion  layers  to  reduce  flooding  [5],  etc.  However, 
there  is  a  still  long  way  to  go  before  solving  such  problems 
satisfactorily,  At  the  same  time,  researchers  have  turned  to  devel¬ 
oping  PEMFCs  operated  at  higher  temperature  (above  120  °C). 
Different  high  temperature  PEMFC  systems  have  been  proposed 
[6-9],  as  reviewed  recently  [10].  One  system  is  the  PBI/H3PO4 
(phosphoric  acid-doped  polybenzimidazole)  high  temperature 
PEMFC.  Researchers  have  focused  their  research  work  on  its 
proton  conductivity  mechanism,  catalyst  activity  and  system 
integration  [11-13].  However,  modeling  of  the  PBI/H3PO4  high 
temperature  PEMFC  is  not  well  documented,  Recently,  Ched- 
die  et  al.  [14,15]  reported  a  one  dimensional  mathematic  model 
of  the  high  temperature  PBI/H3PO4  system,  but  their  study 
was  limited.  In  this  study,  a  two  dimensional  model  based  on 
electrochemical  methods  has  been  constructed,  the  LSY  and 
EIS  electrochemical  techniques  supplied  to  the  model  with  the 
cathode  exchange  current  density  and  cell  internal  resistance, 
generate  good  agreement  between  simulated  and  experimental 
results.  Our  model  offers  a  new  way  to  predict  the  performance 
of  a  PBI/H3PO4  high  temperature  PEMFC  unit  cell  based  on 
electrochemical  techniques. 
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Nomenclature 

A  MEA  active  area 

C  concentration 

d  thickness 

E  voltage 

F  Faraday  constant 

i  current  density 

K  gas  permeability 

N  mass  flux 

Q  charge 

R  universal  gas  constant 

S  entropy 

T  temperature 

u  velocity 

Greek  letters 
a  B-V  coefficient 

s  gas  porosity 

y  B-V  coefficient 

0  overpotential 

H  gas  viscosity 


2.  Experimental 

2.7.  Electrode  and  MEA  preparation 

Catalyst  mixtures  were  prepared  by  directly  mixing  TKK 
46.2%  Pt/C  catalyst-cell  A  (TKK  29%  Pt/C  catalyst-cell  B  and 
J-M  20%  Pt/C  catalyst-cell  C),  NMP  solution  of  PBI  and  PVDF. 
The  mixtures  were  stirred  thoroughly  before  they  were  applied 
to  the  wet-proofed  carbon  paper  (SGF  GDF  31  BC,  300  pan 
of  thickness)  by  a  tape-casting  technique.  The  electrode  fab¬ 
rication  was  completed  by  drying  at  80  °C  for  2h  and  130  °C 
for  0.5  h  in  air.  The  platinum  loading  in  the  electrode  of  cells 
A-C  was  about  0.8  mg  cm-2,  0.75  mg  cm-2  and  0.78  mg  cm-2, 


respectively.  The  same  electrode  was  used  both  for  the  anode 
and  the  cathode.  The  PBI  membrane  (35  ptm,  provided  by  our 
international  research  partner)  was  first  doped  with  85%  H3PO4 
at  60  °C  to  obtain  proton  conduction.  The  unit  cell  was  prepared 
as  follows:  firstly,  the  doped  PBI  membrane  was  sandwiched  in 
between  two  electrodes,  and  then  placed  between  two  graphite 
plates  with  parallel  gas  channels.  Two  aluminum  end  plates  with 
attached  heaters  were  used  to  clamp  the  graphite  plates  by  eight 
bolts  with  a  torque  of  3  N  m.  The  active  single  cell  area  was 
7.8  cm2. 

2.2.  Polarization  curves  measurement 

The  unit  cell  was  installed  on  a  fuel  cell  test  platform.  The  cell 
was  operated  at  150  °C,  0.2  MPa  pressure  (absolute  pressure) 
was  applied  at  both  anode  side  and  cathode  side,  and  the  flow 
rate  was  100  mFs-1  over  both  sides.  Since  the  gases  were  not 
humidified  before  entering  the  cell,  the  water  content  in  the  gas 
stream  was  negligible.  The  voltage  was  recorded  when  it  reached 
steady  state  at  different  currents. 


Fig.  2.  Hydrogen  desorption  oxidation  curves  for  unit  cells. 
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Fig.  1.  (a)  Equivalent  circuit  and  (b)  experimental  (dots)  and  simulated  (line)  AC  impedance  diagrams  for  unit  cells. 
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Gas  Diffueion  Layer 


Fig.  3.  The  computational  region.  Half  flow  channel  =  0.5  mm;  width  of  shoulder  =  1  mm;  GDL  thickness  =  300  [xm. 


2.3.  AC  impedance  measurement 

The  AC  impedance  measurement  was  conducted  on  the  unit 
cell  operated  at  constant  load  (600  mA  cm-2).  The  test  fre¬ 
quency  range  was  from  1  kHz  to  100  mHz  with  logarithmic 
spacing  and  the  AC  amplitude  was  10  mV.  The  electrochemi¬ 
cal  work  station  was  an  EG&G  model  1025  FRA  and  model 
263 A.  The  equivalent  circuit,  measured  AC  impedance  diagrams 
and  equivalent  circuit  simulations  are  shown  in  Fig.  l(a  and  b). 
In  the  equivalent  circuit  shown  in  Fig.  1(a),  L  represents  the 
inductance  element;  R^f  represents  high  frequency  resistance 
that  reflects  the  internal  resistance  of  the  unit  cell;  Rct  represents 
charge  transfer  resistance  that  reflects  the  resistance  of  chemical 
reaction,  it  is  widely  accepted  that  the  anodic  reaction  resistance 
is  negligible,  so  Rct  stands  for  the  charge  transfer  resistance  of 
cathodic  reaction;  CPE  is  the  constant  phase  element  that  used  to 
simulate  the  porous  surface  of  the  electrode.  Fig.  1(b)  shows  that 
the  simulated  AC  impedance  results  agree  with  the  experimental 
data  very  well,  so  the  proposed  equivalent  circuit  in  Fig.  1(a)  is 
reasonable. 

2.4.  Unit  cell  LSV  test 

The  FSV  test  was  measured  in  a  three-electrode  configura¬ 
tion:  WE  was  cathode  side  (purged  by  N2  with  flow  rate  of 
100  mFs-1),  anode  side  acted  as  RE  (DHE)  and  CE  simulta¬ 
neously  (purged  by  H2  with  flow  rate  of  100  mF  s-1),  the  FSV 
curves  were  recorded  from  0  V  (versus  DHE)  to  0.6  V  (versus 
DHE)  with  the  scan  rate  of  20  mV  s-1.  Typical  unit  cell  FSV 
curves  are  listed  in  Fig.  2. 

3.  Model  description 

3.1.  Model  assumptions 

The  following  assumptions  are  made  in  the  development  of 
the  model: 

(1)  Steady  state  operation. 

(2)  Even  temperature  distribute  in  electrode. 

(3)  The  reactant  gases  are  ideal  gases. 

(4)  Isotropic  macro-homogeneous  porous  diffusion  phase. 


(5)  Anode  overpotential  is  omitted. 

(6)  Catalyst  treated  as  interface  rather  than  a  region. 

3.2.  Computational  region  description 

The  computational  region  was  showed  in  Fig.  3,  as  in  Fig.  3, 
the  cross  region  below  two  half  flow  channels  and  one  shoulder 
was  selected. 

3.3.  Brief  introduction  to  govern  equations 
Overall  voltage 

E  =  Ev  77  ohm  0  c  (1) 


Table  1 

Physical  parameters  and  basic  conditions  [14-17] 


Items 

Value 

ME  A  active  area  A  (cm1 2) 

7.8 

EAS — electrochemical  area  surface  (cm2  mg- 1  Pt) 

As  calculated 

Active  area  factor,  Farea 

As  calculated 

Flow  channel  half  width,  a  (m) 

5  X  10“4 

Shoulder  width,  b  (m) 

10"3 

Oxygen  molar  concentration  at  flow  channel, 

56.8 

Co2,o  (mol  m-3) 

Water  molar  concentration  at  flow  channel, 

0 

Ch2o,o  (molm-3) 

Oxygen  reference  concentration,  Co2,ref  (molm-3) 

3.39 

GDL  thickness,  d\  (m) 

3  X  10“4 

Catalyst  layer  thickness,  J2  (m) 

5  x  10“5 

O2  effective  diffusion  coefficient,  Dq3  (m2  s-1) 

1.2  x  10“6 

Water  effective  diffusion  coefficient,  D^0  (m2  s-1) 

1.4  x  10“6 

Cathode  permeability-coefficient,  kp?c  ( m 2) 

5.5  x  10“15 

ORR  exchange  current  at  planar  Pt  surface, 

1(T3 

if), area-planar  (A  m  ) 

Volume  exchange  current,  z'o.voiume  (Am-3) 

As  calculated 

Working  pressure,  P  (Pa) 

2  x  105 

High  frequency  resistance  (internal  resistance),  Rm  (£2  cm-2) 

As  measured 

Working  temperature,  T  (K) 

423 

Cathode  porosity,  sc 

0.4 

Pt  loading,  Lpt  (mg  cm-2) 

1 

Cathode  gas  viscosity,  /ic  (kgm-1  s-1) 

2  X  1(T5 

Hydrogen  atom  absorption  factor:  Qre f  (|xQcm-2) 

220 

a  a 

0.5 

«c 

1 

Yc 

1 
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Oohm  —  I  ’  ^hf 


where  R^f  can  be  determined  by  EIS 


2  RT 

T0)  H - —  •  In 

nF 


(2)  The  mass  flux  of  oxygen  and  gaseous  water  can  be  divided 

into  two  terms  as: 


Nq2  —  — £:jDq^VCo2  +  sCq2u 


^  Nh2o  =  -£Dft20^Cn2o  +  £Cu2ou 


where  reversible  cell  potential  is  1.185  at  T=7o  =  298K  and 
reactant  gases  at  its  standard  conditions. 

The  following  relation  is  valid  for  400  K  <  T<  500  K: 

—  =  -18.449  -  0.01283  •  T  (4) 

n 


The  first  term  —  £/)effVC  in  Eqs.  (7)  and  (8)  represents  diffu¬ 
sion  flux  and  the  second  term  sCu  in  Eq.  (7)  and  (8)  represents 
convection  flux. 

If  combined  Darcy  law-equation  (9)  which  is  typically  used 
to  descript  porous  diffusion  media  into  Eqs.  (7)  and  (8),  Eqs.  (5) 
and  (6)  can  be  written  in  the  form  as  Eqs.  (10)  and  (11) 


In  cathode  diffusion  layer 

V^Vo2  =  0 
V%o  =  0 


(5) 


(6) 


kvRT 

U  =  - V(Co2  +Ch2o) 

/J 


(9) 


,  pff  £Co2kvRT 

V  (  -£D^  VCo2 - - V(Co2  +  Ch2o)  )  =  o  (10) 


Given  i]c  — cathode  oveipotential 


▼ 

Calculate  02  anclH20  concentration  distribution  according  to 

equation  (10)  and  (11) 


▼ 

Calculate  ciurent  distribution  according  to  equation  (16)  and 
Current  density  according  to  equation  (17) 


T 

Calculate  ciurent  denisty  according  to  equation  (1) — (4) 


Change  cathode  oveipotential  and  repeat  above  procedure 
The  polarization  curves  were  simulated 


Fig.  4.  Solution  algorithm. 
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(  eff  sC^okryRT  \ 

V  ( -f<oVCH2o - — V(Co2  +  CH2o)  )  =  0 

(11) 

3.4.  Boundary  conditions 

Boundaries  2  and  3  are  symmetric  interface,  so  the  mass  flux 
at  normal  direction  of  O2  and  H2O  is  zero: 

No2,n  =  0;  NU2o,n  =  0  (12) 

Boundaries  4  and  5  are  the  interface  between  flow  channel 
and  GDL,  so  the  constant  concentration  is  applied  as: 

Cq2  =  Cq2,o;  Gh2o  =  Gh2o,o  (13) 


Table  2 


Cathode  volume  exchange  current  density  for  unit  cells 


Cell  number 

Cell  A 

Cell  B 

Cell  C 

Cathode  volume  exchange 

1649 

1404.4 

996 

current  density  (Am-3) 

Table  3 

Internal  resistance  for  unit  cells 


Cell  number 

Cell  A  Cell  B 

Cell  C 

Internal  resistance  (£2  cm2) 

0.209  0.218 

0.21 

Boundary  6  is  the  interface  below  shoulder,  the  mass  flux  at 
normal  direction  of  O2  and  H2O  is  zero: 

No2,n  =  0;  ^VH2o,n  =  0  (14) 

Boundary  1  is  the  catalyst  layer  interface,  where  oxygen  is 
consumed  and  water  is  formed,  so  the  boundary  conditions  can 


be  described  as  below  combined  with  the  general  accepted  elec¬ 
trochemical  B-V  Eq.  (16): 


Nu20,n  = 


I' 

IF 


(15) 


CO 

1 

57 

E 

• 

0  47.5 

E 

38 

O 

2 

28.5 

c 

0 

O 

19 

c 

0 

0 

9.5 

CM 

0 

O 

0  57 

1 

E 

o  45.6 
.1  34.2 

CO 

22.8 

0 

o 

8  11-4 

o 


CM 

o 


0 


Fig.  5.  Oxygen  concentration  surface  plots  at  different  cathode  overpotentials  for  cell  A:  (a)  ijc  =  0.05  V;  (b)  ijc  =  0.15  V;  (c)  rjc  =  0.25  V;  (d)  ijc  =  0.35  V. 
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where 


/'  =  Hi 


Co¬ 


re 


area 


^02,  ref 


exp 


oicFif 

RT 


exp 


<XaF% 


RT 


(16) 


3.5.  Numerical  values 


The  operating  conditions  and  electrode  parameters  are  shown 
in  Table  1 . 

3.6.  Solution  algorithm 

Fig.  4  shows  the  solution  algorithm  used.  After  the  current 
density  along  the  catalyst  layer  is  obtained,  the  current  density 
can  be  calculated  by  the  following: 


*L 


I  =  L 


r  dx 


(17) 


According  to  the  solution  algorithm  showed  in  Fig.  4,  the  polar¬ 
ization  curves  can  be  simulated.  The  equations  of  this  model 
were  solved  using  commercially  available  program  FEMLAB® . 

3.7.  Parameters  derived  from  electrochemical 
measurements 

The  cathode  volume  exchange  current  density  was  calculated 
according  to  the  following: 


EAS  = 


Qh 


Gref  ’  ^pt  •  A 


-^area  —  EAS  •  Lpt 

h),area  =  EAS  •  Lpt  •  z’o, area-planar 

*0,  volume  =  EAS  •  Lpt  •  4), area-planar  '  ^2 


(18) 

(19) 

(20) 
(21) 


Eq.  (19)  gives  the  active  area  factor  that  represents  the  ratio 
between  electrochemical  area  surface  and  apparent  electrode 
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Fig.  6.  Gaseous  water  concentration  surface  plots  at  different  cathode  overpotentials  for  cell  A:  (a)  i;c  =0.05  V;  (b)  i;c  =0.15  V;  (c)  i;c  =  0.25  V;  (d)ijc  =  035  V. 
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area;  Eqs.  (20)  and  (21)  give  the  cathode  area  exchange  current 
density  and  volume  exchange  current  density  respectively.  The 
calculated  cathode  volume  exchange  current  densities  for  differ¬ 
ent  unit  cells  are  shown  in  Table  2.  The  electrochemical  method 
offers  a  new  way  to  determine  the  cathode  area  exchange  current 
density  and  volume  exchange  current  density,  the  value  of  cal¬ 
culated  cathode  volume  exchange  current  density  is  comparable 
to  that  generally  cited  in  literature  [17] — 440  Am-3. 

As  to  cell  internal  resistance,  the  value  was  determined  by  the 
AC  impedance  measurement  based  on  equivalent  circuit  anal¬ 
ysis,  no  obvious  mass  transfer  resistance  related  to  the  transfer 
phenomenon  caused  by  liquid  water  [18]  in  the  lower  frequency 
region  was  found  and  no  liquid  water  exists  in  the  high  temper¬ 
ature  PEMFCs,  only  the  internal  resistance  R^f  and  the  charge 
transfer  resistance  Rci  appear  in  the  AC  impedance  diagram,  so 
it  is  reasonable  to  assume  that  ohmic  polarization  and  cathode 
reaction  polarization  are  the  main  contributions  to  the  polariza¬ 
tion,  as  is  shown  in  Eq.  (1). 

The  Rm  for  different  cells  is  listed  in  Table  3.  It  should  be 
noted  that  the  R^f  in  Table  3  is  the  internal  resistance  when  the 
cell  is  operated  at  the  current  density  of  600  mA  cm-2,  in  gen¬ 
eral,  the  Rtf  dependence  on  current  density  should  be  measured 
for  Rtf  is  dependant  on  current  density,  however,  in  our  exper¬ 
iment,  we  found  that  the  internal  resistance  derived  from  AC 


(a)  X  Position  at  Catalyst  Layer  (y=-0.00015  m)/m 


impedance  is  almost  stable  when  changing  the  load,  we  attribute 
this  phenomenon  to  the  proton  conduction  of  PBI/H3PO4  mem¬ 
brane  which  is  independent  of  the  water  content.  For  this  reason, 
the  internal  resistance  of  the  unit  cell  is  taken  as  a  constant. 

4.  Simulation  results 

Fig.  5(a-d)  list  the  O2  concentration  distributions  at  different 
cathode  overpotentials  for  cell  A  (cell  B  and  cell  C  have  similar 
results  that  is  omitted  here)  in  GDL  based  on  the  two  dimensional 
model  described  above,  from  Fig.  5,  the  O2  concentration  of  the 
region  near  catalyst  layer  is  lower  than  that  near  flow  channel 
for  O2  is  consumed  at  catalyst  layer;  and  the  O2  concentration 
gradient  increases  with  increase  of  cathode  overpotential,  the 
mean  O2  concentration  at  the  catalyst  layer  decreases  to  21.9 
when  the  cathode  overpotential  is  0.35  V  (Fig.  5(d));  in  addi¬ 
tion,  the  O2  concentration  of  the  GDF  region  under  shoulder 
(—0.5  mm  <  v  >  0.5  mm)  is  lower  than  that  of  the  GDF  region 
under  the  flow  channel  for  the  GDF  region  since  O2  can  only 
enter  this  region  from  the  X  direction. 

Fig.  6(a-d)  list  the  gaseous  water  concentration  distribu¬ 
tions  at  different  cathode  overpotentials  for  cell  A  based  on 
the  two  dimensional  model  (cells  B  and  C  have  similar  results 
but  are  omitted  here),  from  Fig.  6,  the  gaseous  water  concen- 


(k)  X  Position  at  Catalyst  Layer  (y=-0.00015  m)/m 


(d)  X  Position  at  Catalyst  Layer  (y=-0.00015  m)/m 


Fig.  7.  Current  density  distributions  at  catalyst  interface  at  different  cathode  overpotentials  for  cell  A:  (a)  r)c  =  0.05  V;  (b)  ijc  =  0.15  V;  (c)  rjc  =  0.25  V;  (d)  r)c  =  0.35  V. 
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Table  4 

Current  density  under  shoulder  and  under  flow  channel  at  different  cathode  overpotential 
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Cathode  overpotential 

0.05 

0.15 

0.25 

0.35 

Mean  current  density  under  shoulder  (A  cm  2) 

0.0513 

0.183 

0.522 

0.960 

Mean  current  density  under  flow  channel  (A  cm-2) 

0.0521 

0.195 

0.689 

1.560 

tration  of  the  region  near  catalyst  layer  is  higher  than  that  of 
the  region  near  flow  channel  for  water  that  is  produced  by 
the  catalyst  layer;  and  the  water  concentration  gradient  also 
increases  with  increase  of  cathode  overpotential.  In  addition, 
the  water  concentration  of  the  GDL  region  under  the  shoulder 
is  higher  than  that  of  the  GDL  region  under  the  flow  channel 
because  the  transport  of  water  produced  at  the  catalyst  layer  is 
blocked  by  the  shoulder  that  does  not  contact  the  flow  channel 
directly. 

Fig.  7(a-d)  show  the  current  density  distributions  along  the 
catalyst  layer  for  cell  A  at  different  cathode  overpotentials,  also 
the  equivalent  plots  of  cell  B  and  cell  C  are  omitted  in  this  paper 
for  brevity.  The  solid  line  in  Fig.  7  represents  the  Gaussian  fit 


of  the  current  density  distribution  along  the  catalyst  layer,  and 
the  obtained  continuous  fitted  line  is  used  to  calculate  the  mean 
current  density  along  catalyst  layer  according  to  Eq.  (17).  From 
Fig.  7  the  current  density  distributions  along  catalyst  layer  are 
proportional  to  the  O2  concentration  distribution  along  the  cat¬ 
alyst  layer  (please  refer  to  Fig.  5),  and  the  current  density  at  the 
catalyst  layer  under  the  shoulder  is  lower  than  that  under  the 
flow  channel,  this  effect  is  more  obvious  when  the  cathode  over¬ 
potential  increases  (Table  4).  For  this  reason,  the  catalyst  under 
the  shoulder  is  not  fully  utilized  especially  when  the  cathode 
overpotential  is  high. 

Fig.  8(a-c)  show  the  experimental  and  simulated  polarization 
curves  in  this  study,  from  Fig.  8,  the  experimental  and  simulated 


Fig.  8.  Experimental  and  simulated  polarization  curves:  (a)  cell  A;  (b)  cell  B;  (c)  cell  C. 
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results  agree  with  each  other  very  well,  which  verifies  the  relia¬ 
bility  of  our  model. 

It  should  be  noted  that  this  model  is  a  preliminary  two  dimen¬ 
sional  model  study  of  the  PBI/H3PO4  high  temperature  PEMFC 
fuel  cell,  and  the  catalyst  layer  is  treated  as  a  reaction  interface, 
so  it  cannot  reflect  the  effect  of  catalyst  microstructure  on  the 
cell  performance,  more  detailed  model  research  work  combined 
with  electrochemical  techniques  is  under  development  in  our 
group. 

5.  Conclusion 

A  two  dimensional  cathode  model  combined  with  elec¬ 
trochemical  data  was  constructed  in  this  study  based  on  the 
PBI/H3PO4  high  temperature  PEMFC  fuel  cell,  the  oxygen  and 
gaseous  water  concentration  distributions  were  calculated  based 
on  the  model,  the  polarization  curves  of  three  unit  cells  were  sim¬ 
ulated,  and  the  simulated  and  experimental  results  agree  well. 
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